Quantum phase transitions in the Bose-Fermi Kondo model 
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We study quantum phase transitions in the Bose-Fermi Kondo problem, where a local spin is coupled 
to independent bosonic and fermionic degrees of freedom. Applying a second order expansion in 
the anomalous dimension of the Bose field we analyze the various non-trivial fixed points of this 
model. We show that anisotropy in the couplings is relevant at the SU(2) invariant non Fermi liquid 
fixed points studied earlier and thus the quantum phase transition is usually governed by XY or 
Ising-type fixed points. We furthermore derive an exact result that relates the anomalous exponent 
of the Bose field to that of the susceptibility at any finite coupling fixed point. Implications on the 
dynamical mean field approach to locally quantum critical phase transitions are also discussed. 
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I. INTRODUCTION 

Quantum phase transitions in quantum-impurity mod- 
els such as the two impurity Kondo modeljd □ anisp,taopic 
Kondo modelsfi multi-channel crystal field modelsBQ and 
the Bose-Fermi Kondo model (BFKM) attracted recently 
a lot of interest. These models provide the simplest 
realizations of zero temperature 'quantum phase tran- 
sitions': They have several possible strongly correlated 
ground states, each described by some zero temperature 
(T = 0) critical point, and display a phase transition be- 
tween these as one changes the coupling constants that 
control their behavior. The quantum critical point of 
these systems is in many cases described by an unstable 
non-Fermi liquid impurity model with singular thermo- 
dynamic and transport properties. Recently, it became 
also possible to study some of these quantum-phase tran- 
sitions experimentally in mesoscopic deviceso. 

The present paper is devoted to a detailed analysis of 
the Bose-Fermi Kondo problem, consisting of a local im- 
purity spin, S a , that couples both to fermionic (ip) and 
bosonic (<f>) fields. In the present paper we restrict our 
discussion to the case of a spin 5 = 1/2 local moment, 
but our results can be trivially generalized for larger spins 
and many electron channels. The bosonic degrees of free- 
dom usually represent a fluctuating magnetic order pa- 
rameter field at a quantum critical point p 113 and exhibit 
power-law correlations-at zero temperature, while ip rep- 
resents free Fermions.li-3 This model has been proposed 
in Ref. [l5| to describe the 'locai-nuantum phase transi- 
tion' in alloys like CeCu5_ x Au x £j lill. Its purely bosonic 
version emerges in the context of spin-glass mean field 
theories,Bli3 and it has been proposed in Ref. [l3| to de- 
scribe non-magnetic impurities in two-dimensional quan- 
tum antiferromagnets. 

The Bose-Fermi Kondo problem in its full version has 
been analyzed independently by Smith and Sijij and 
Senguptalij by performing a leading order expansion in 
the anomalous dimension e of the Bose field, <f>. It has 
been shown in Refs. [ll],[l2| that the competition between 
the bosonic and fermionic fields gives rise to a quantum 



phase transition, with the two stable quantum phases 
corresponding to purely bosonic and purely fermionic 
models, respectively. At the 'fermionic' fixed point the 
impurity spin is screened by a Kondo effect, and the 
bosonic field decouples. At the bosonic fixed point, on 
the other hand, fermionic fields are irrelevant, the impu- 
rity spin becomes partially screened and the properties of 
this fixed point are coiitrolled by the anomalous dimen- 
sion e of the Bose field.Ej An expansion to second order in 
e has been performed for the isotropic and purely bosonic 
model in Ref. |l^. 

In the present paper, we study the fully anisotropic 
Bose-Fermi Kondo model, described by the interaction 
Hamiltonian: 



Hint = H 



K 



(1) 



Here A a and j a (a — x, y, z) denote dimensionless cou- 
pling constants, and A is a high energy cutoff. We as- 
sume that the fermionic imaginary time propagator cor- 
responds to a free degenerate Fermi gas and it therefore 
decays as ~ 1/r at T = 0: 



2tt 



(2) 



while the bosonic field shows critical imaginary time cor- 
relations at zero temperature with an anomalous dimen- 
sion e 



(T0 Q (r)^(O)) = est. 







(3) 



In this work we consider e as a given external parameter 
which will always be taken to be positive. In practice, it 
can be generated by the critical dynamics of some, spin 
degrees of freedom at a quantum critical point ,E2fli3 and 
for a dense system of magnetic impurities it must be de- 
termined self- consistently.EJ 
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To analyze the quantum phase transition we carry out 
a second order expansion in e. As already remarked in 
Rcf. [12], anisotropy in the couplings A" and 7" is rele- 
vant, even in the purely bosonic model: (a) We find that 
the stable bosonic phase is usually anisotropic, and (b) 
the quantum phase transition in the Bose-Fermi Kondo 
model is in many cases governed by the new anisotropic 
fixed points discussed below and not the SU(2) invari- 
ant fixed point discussed in detail in Refs. |l^Ju] , [il| We 
also obtain some exact results for the quantum critical 
behavior that follow from a Ward identity and are valid 
to all orders in e. These results parallel the exact results 
of Refs. |l3||l§| obtained for the bosonic versions of the 
model, and raise important questions in the context of 
the theory of Ref. [15]. 

Postponing the detailed and rather involved computa- 
tions to the following sections, let us briefly discuss our 
main results and the methods applied here. 

To describe the quantum phase transitions within the 
Bose-Fermi Kondo model (BFKM) we used the powerful 
machinery of multiplicative renormalization group (RG) . 
The infinitesimal renormalization group transformations 
are most conveniently described by scaling equations, 
which can be obtained by performing a perturbative ex- 
pansion of the various vertex functions in the couplings 



Xf, and 



and then gradually eliminating high- 



energy degrees of freedom by reducing the cutoff A — * A'. 
The reduction of the cutoff must be compensated by re- 
scaling the couplings in such a way that the physical ob- 
servables and the singular part of the free energy remain 
invariant. These calculations, detailed in the main body 
of the paper, lead to the following Gell-Mann-Low equa- 
tions: 



^=/3i b) ({A 7 },{3 7 }), 



(4) 



where I = hi(4f) denotes the scaling variable with Ao 
the initial value of the cut-off, fla^ and (3a^ stand for the 
fcrmionic and bosonic beta functions, respectively, and 
we introduced new 'gauge invariant' couplings: 

9v = ll ■ 

In fact, since the theory is invariant under the Z2 trans- 
formation 4> ~ > —<fi, physical quantities may only depend 
on ll- 

We computed to 0(e 3 ) by performing a per- 

turbative expansion of the various vertex functions in 
A7 ~ 5j ~ £■ The explicit expressions are given by 
Eqs. ( pi[ ) and (|65|), here we only give the leading order 
results: 



= X y X z - ^X x (g y +g z ) 

dg x , \ 

-jT = e 9x - 9x{g y +gz) ■ 



(5) 
(6) 



The other three equations can be obtained by cyclic per- 
mutations. 

The first term in Eq. (||) corresponds to the usual dy- 
namic screening present in the fermionic Kondo model, 
while the linear term in Eq. (^) reflects the anomalous di- 
mension of the bosonic field <j>. The last terms in Eqs. 
and (^|) are due to the dissipative coupling between the 
impurity spin and the bosons: The bosonic heat bath 
adjusts itself to a given spin configuration and strongly 
suppresses processes with spin flips. 

Eqs. (Q) turn out to have a number of non-trivial small 
coupling fixed points with couplings X^,g^ ~ e: These 
control the possible quantum phases and the transitions 
between them in the BFKM. 



A. Non-trivial fixed points 

Before we turn to the general case, it is instructive to 
discuss the two special cases of (1) purely bosonic cou- 
plings and (2) full SU (2) symmetry with = g and 
A ; , A. 



1. The purely bosonic model 

For A, = the RG flows are shown in Fig. [j] Two 
types of non-trivial fixed points appear: 

(a) Bosonic SU(2) fixed point with g x = g y — g z : This 
is the fixed point analyzed in Ref. [l3|, and suggested 
to control the behavior of non-magnetic impurities in 
a two-dimensional antiferromagpet at the antiferromag- 
netic quantum phase transition.E3 

(b) Three Bosonic XY fixed points of the type g x — g y ^ 
9z = 0. 

Both fixed points are unstable against spin anisotropy. 
This instability is cha ra c te ristic to quantum states with 
finite residual entropy t3t3 The system tries to get rid of 
this entropy by breaking the symmetry, and in fact, the 
only stable fixed points of the purely bosonic model seem 
to be the three infinite coupling Ising fixed points of the 
form g x = g y = ^ g z — ► 00. 
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FIG. 1. Fig. (a) Sketch of the renormalization flows in the 
purely bosonic model. We find one non-trivial SU(2)-invariant 
fixed point and three fixed points with XY symmetry. Both 
are unstable against breaking the SU(2) and XY symme- 
tries, respectively. Fig. (b) Quantum phase transition in the 
purely bosonic model with XY symmetry as the anisotropy, 
$9 = 9z — ff-L changes sign. The dashed lines indicate 
cross-over regions corresponding to 

Even the behavior of the purely bosonic model is al- 
ready very rich. In a realistic number of cross- 
overs may occur between SU(2)-type, XY or Ising behav- 
iors. In particular, for a system with XY-type symmetry 
a quantum phase transition occurs between the XY-type 
and Ising fixed points as Sg = g z — g x — g z — g y changes 
sign (see Fig |l|). In this case, the quantum phase transi- 
tion is controlled by the SU(2) type fixed point below an 
energy scale T*. For small g± = g x = g y as g z = g one 
can determine T* by integrating the scaling equations 
Eq. (^), and is approximately given by 

T* ~ Aog 1 ^. 

As T — > another cross-over occurs at an energy scale 

T** ~ (Sg) Vr T* , 

where y r = e/2+e 2 /2 + ... is the scaling dimension of spin- 
anisotropy at the SU(2) fixed point. Between T* and T** 
the SU(2)-invariant Bosonic model describes the behavior 
of the impurity, while below T** the physical properties 
of the model are controlled by the Ising (5g > 0) or the 
XY-type (Sg < 0) bosonic fixed points, respectively. 

If g x ^ g y then the behavior is even more complicated, 
and eventually two cross-overs may occur in series, cor- 
responding to the SU (2) — > XY and XY — > Ising tran- 
sitions. The cross-over regions are governed in this case 
by the SU(2) and XY-type fixed points, respectively. 

Although anisotropy is relevant at both 0(e) fixed 
points above, in many cases the SU(2)-invariant or XY- 
type fixed points can be-also of physical interest: In the 
case of couprates, e.g.,t!j spin-orbit interaction is weak 
and therefore the anisotropy is presumably weak. In this 
case the SU(2)-symmetrical fixed point may appropri- 
ately describe the physics over a wide energy range of 
interest. On the other hand, in most magnetic heavy 
Fermion compounds spin-spin interactions are generically 
anisotropic ,t£l and the physics is controlled by XY type 
or Ising fixed points. 



2. The SU(2) symmetric Bose-Fermi Kondo Model 

For the sake of clarity, we first discuss the renormal- 
ization group flows in the SU(2) symmetric case with 
9x = 9y = gz = g and X x = \ y = X z = A. The results 
we obtain to second order in e are qualitatively the same 
as the leading order ones derived in Refs. |ll|Jl2| . The 
flows are sketched in Fig. ||. Two stable quantum phases 
appear: 

(a) The Kondo. state (X — > oo, g = 0): This is the familiar 
Kondo statc.EII The local moment is completely screened 
by the Fermions and forms a Fermi liquid. It is thus fully 
decoupled from the bosonic field. At the Kondo fixed 
point both anisotropy and the bosonic field are irrelevant. 

(b) Purely bosonic SU(2) invariant fixed point (g s» e/2, 
X = 0): This is the same fixed point as the one we dis- 
cussed in the previous subsection. At this fixed point the 
coupling to the Fermions is irrelevant, however, break- 
ing the spin anisotropy, as discussed earlier, is a relevant 
perturbation. It is though a stable fixed point if for some 
reason exact SU(2) symmetry is guaranteed, however, it 
is generally unstable in non-cubic systems with spin-orbit 
interaction. At_this fixed point the impurity is only par- 
tially screened £3 

These two phases are separated by an interesting quan- 
tum critical point: 

(c) The SU(2) symmetrical Bose-Fermi fixed point (<? ~ 
A ~ e) . This fixed point governs the quantum phase tran- 
sition between the bosonic and Fermionic SU(2) invariant 
phases in case of full SU(2) symmetry. However, similar 
to the SU(2) invariant Bose fixed point, anisotropy is 
relevant at this fixed point (in fact, it is the leading rele- 
vant operator, i.e., it is more relevant than the operator 
corresponding to the quantum phase transition.) There- 
fore, in general, it is not this fixed point that controls the 
quantum phase transition. 




SU(2) Bose 6 I 

FIG. 2. Renormalization group flows for the SU(2) sym- 
metrical Bose-Fermi Kondo model. "BF" denotes the SU(2) 
symmetrical Bose-Fermi fixed point. Here the spin is partially 
screened and both the Fermi and Bose fields couple to it. At 
the Bose fixed point the fermionic degrees of freedom fully 
decouple from the spin. 
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3. The general case: Anisotropic BFKM 

The general anisotropic Bose-Fermi Kondo model has 
all the fixed points discussed above. In addition, however, 
we find two new types of fixed points with XY and Ising 
symmetry: 

(a) The XY Bose-Fermi fixed point of the type \ x — 
A y 7^ A 2 and g x — 9y ^ 9z = 0: Similar to the SU(2) 
Bose-Fermi fixed point, this fixed point separates the XY 
or Ising bosonic phases and the SU(2) invariant Kondo 
phase, and controls the quantum phase transition be- 
tween them in case of XY symmetry. 

Breaking the SU(2) symmetry is a relevant perturba- 
tion at the XY-Bose-Fermi fixed points too: It corre- 
sponds to an operator within the critical surface sepa- 
rating the bosonic and fermionic quantum phases. If we 
break the XY symmetry in this direction, the RG flows 
stay within the critical surface and end up in a new strong 
coupling fixed point: 

(b) The Ising Bose-Fermi fixed point with structure, 
Xx — X y ^ Az and g x = g y = ^ g z , with g z , \ z ~ 1. 

This fixed point is outside the region of applicabil- 
ity of our perturbative RG. We believe, however, that 
this strong coupling fixed point is not a mere artifact 
of our approximations, and that this is the fixed point 
that should describe the phase transition in case of Ising 
symmetry. 

The complicated structure of the flows and the con- 
nection between the various fixed points is summarized 
in Fig. | 
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FIG. 3. Summary of various fixed points and the possible 
flow directions. 



B. Physical quantities 

Within the e-expansion, we can compute various phys- 
ical quantities. The scaling properties of these are de- 
termined by the scaling dimensions and structure of the 
various relevant and irrelevant operators. Instead of enu- 
merating all of these at each fixed point we rather list 
only the most important scaling dimensions in Table 0. 



At the 0(e) fixed points we enumerated three scaling 
dimensions: 

(a) The leading relevant operator. This is usually associ- 
ated with breaking the symmetry of the fixed point. 

(b) We also gave the scaling dimension of the most rele- 
vant operator respecting the symmetry of the fixed point. 
This operator at the Bose-Fermi fixed points corresponds 
to the quantum phase transition between the bosonic and 
fermionic fixed points with that symmetry. 

(c) Finally, we also enumerated the scaling dimension of 
leading irrelevant operators, which determine the scaling 
properties of physical quantities at the quantum criti- 
cal point (or the symmetry-stabilized XY or SU(2) Bose 
fixed points). 

1. Susceptibility 

Magnetic field is a relevant perturbation in all non- 
trivial 0(e) fixed points. To study its effect we added 
the following term to the Hamiltonian: 

H m = — huS a , (7) 

a— x : .. : z 

where h a denote the components of a local magnetic field. 
As we shall prove later, the scaling equations of the di- 
mensionless magnetic field, h a = h a /A, are determined 
by the same beta- functions as those of the vertex g a : 

^T = l-| + ^ ] «M. i9,}) + 0(h) . (8) 

This remarkable result, which is very similar to the one 
obtained in Ref. [l8| and follows from a Ward identity, has 
important consequences: By Eq. (^) it implies that at a 
fixed point with coupling g* a ^ 0, the scaling dimension of 
the a component of the local magnetic field only depends 
on the anomalous dimension of the bosonic bath, and is 
exactly given by yn.a = 1 — e/2. 

As explained in Ref. |2l],e.<7., the dimension xq of a 
scaling operator at a non-trivial fixed point is defined 
through the asymptotic behavior of its correlation func- 
tion: (0(r)0(O)) ~ l/r 2x °, where r denotes the coor- 
dinate variable. The dimension xq is directly related 
to the scaling dimension yn of the dimensionlcss field 
H that couples to the the scaling operator: Since in a 
D-dimensional theory the free energy has a dimension 
D, yn = D — xq. In our case the impurity spin lives 
in an effective D = + 1-dimensional space, therefore 
this relation immediately relates the scaling dimension 
of the magnetic field h to that of the impurity spin S as 
xs = 1 — Vh- As a consequence, the spin correlation func- 
tion at the SU(2)-invariant Bose and Bose-Fermi fixed 
points decays at T = temperature as 

xl U {2) (r) = (TS a (r)S (O)) = S a0X (r) , (9) 

X(r) = (10) 
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with T* , a dynamically generated energy scale similar to 
the Kondo scale. This result has been derived earlier for 
the SU(2)-invariant bosonic fixed point in Ref. [ p"3| . 

At the XY fixed points, however, the scaling properties 
of the magnetic field depend on its direction. To be spe- 
cific, let us assume that the value of the couplings at the 
fixed point is g* — g* ^ g* — 0. In this case, spin corre- 
lations of the x and y spin components decay at T = 
as 



x\ (t) 



1 



(11) 



while the decay of correlations in the z-direction is fixed 
point specific: 



X 



(XY) 



(r) 



1 



(12) 



where, up to second order in e, r\ z is given by i] z = 2e + e 2 
and rj z = 2e + 7e 2 /4 for the XY Bose and XY Bose-Fermi 
Kondo fixed points, respectively. Note that Eqs. ( fiol) and 
( pd| ) are valid to all orders in e, while the exponents rj z 
can be computed only approximately. 

The results above apply only at T = temperature. 
On can, however, obtain the finite temperature form of 
the correlation functions and their Fourier transforms as- 
suming— conformal invariance at these non-trivial fixed 
points E3 For the imaginary part of the Fourier trans- 
form of the local spin-spin correlation functions, measur- 
able through neutron scattering experiments, we obtain 
the following scaling form: 



Xc 



T 



where the scaling function f(jj,x) behaves as 



C> \x\ 
C< \x\ 



M » T, 
M <T. 



(13) 



(14) 



The function f(r],x) can be expressed analytically from 
the exact expression Eq. (Kq). 



2. Resistivity 

Here we only discuss the impurity resistivity at the 
quantum critical points, i.e., we assume that all relevant 
variables and symmetry breaking terms have been tuned 
to zero or vanish for reasons of symmetry. In this case, 
the temperature dependence of the resistivity is deter- 
mined by the scaling dimension yi < of the leading 
irrelevant operator. We enumerated this latter for the 
various fixed points in Table [j]. 

The simplest way to obtain the resistivity is to com- 
pute the impurity scattering rate. This turns out to be 
proportional to J2 a ^a(T), where A a (T) is the effective 
coupling at energy scale to ~ T that can be computed 



from the scaling equations. Since at the bosonic fixed 
point the effective couplings A Q vanish as ~ A Q (T) ~ T A 
with a power A = — y^ the resistivity there also vanishes 
as 



Rf™(T)~AT 



2A 



(15) 

and the impurity fully decouples from the electrons at 
T = 0. 

The situation at the SU(2) and XY Bose-Fermi fixed 
points is, on the other hand, different. There the cou- 
plings scale to a finite value, corresponding to a finite 
scattering rate, and the leading corrections to this scale 
as 



BZ V {T)~A + BT' 



(16) 



similar to the multichannel Kondo model. 

The above results only apply for the case of a dilute 
system of independent impurities. For the dense sys- 
tems the interaction between individual impurities must 
be taken into account which leads in general to a different 
temperature dependence. 



C. Implications for a dynamical mean field theory of 
locally quantum critical phase transitions 

It has been shown in a beautiful series of neutron scat- 
tering measurements by Schroder et al. that the quantum 
critical behavior in CeCu5_ x Au x is entirely local within 
experimental accuracy. This material shows a quantum 
phase transition from a non-magnetic heavy Fermi liquid 
state to a metallic antiferromagnetic state under Au dop- 
ing. The Fermi liquid energy scale below which the Ce 
spins are screened appears to go to zero as one approaches 
the quantum critical point from the metallic side, indi- 
cating that at the quantum critical point (QGE) spins re- 
main asymptotically free at low temperature. til Schroder 
et al. also find that at the quantum critical point the 
dynamic susceptibility behaves with a good accuracy as 

X-V,T,q) ~ cst.fr-^q) + (T — iuo) a ) (17) 

with an exponent a i=s 0.75, where x~ is an approxi- 
mately frequency and temperature independent function 
with minima at the ordering wave vectors. 

Clearly, the quantum phase transition in CeCu5_ x Au x 
is driven by the competition between two mechanisms 
with which the system can get rid of the residual entropy 
of the Ce impurity spins: The Kondo effect that screens 
magnetic impurities individually, and the magnetic or- 
dering, which tends to align spins and thus destroys the 
Kondo effect. 

Expression ( |l7| ) is consistent with the presence of 
mean-field-like spatial correlations but simultaneous non- 
trivial local correlations in the time direction. Based on 
these observations J13 Si et al. proposed-a self-consistent 
version of the Bose-Fermi Kondo modeLJ as a good can- 
didate to describe the above quantum optical behavior 
within a dynamical mean field approach.!^ 
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1. SU(2) invariant theory 

In their work, Si et al. assumed that the system is 
SU(2) symmetrical in spin space, and therefore identi- 
fied the QCP with the SU(2) invariant Bose-Fermi fixed 
point. In this theory, the exponent e must be determined 
self-consistently, and it is fixed by the condition that the 
correlation functions (TS a (T)S a (0)} and (T0 a (r)0 a (O)) 
must decay asymptotically with the same exponent. 

As we have discussed above, at the SU(2) invariant 
Bose-Fermi fixed point (S(r)S(O)) ~ l/r e as a conse- 
quence -aLa Ward identity. Similar to the case of spin 
glasses,EfE3 this immediately determines the only possi- 
ble exponent: 



SU(2) _ 
e QCP — 1 



(18) 



This result has dramatic consequences: 

(a) It implies that - if a stable self-consistent solution of 
the extended dynamical mean field theory exists - then 
local spin correlations decay as 1/r at the quantum crit- 
ical point and thus the the local susceptibility must di- 
verge logarithmically as a function of frequency or tem- 
perature. This may be consistent with the mean field 
expression ( p"7| ) only if the magnetic fluctuations are ap- 
proximately two-dimensionalBB In this case the exponent 
a would be non-universal and would depend on the spe- 
cific properties of the model. 

(b) As a consequence, the scaling dimension of the local 
magnetic field (i.e. a field applied only to the impurity) 

would be exactly y h — 1/2. This in turn immediately 
implies an h 2 /T scaling, which is clearly in disagreement 
with the experimental data, showing an almost perfect 
h/T scaling. However, in the experiments one applies 
a global field, and we know cases where the the scaling 
dimension of the global field is diffcrcuLitpm that of the 
local field acting only on the impurity.Ejnfl It is therefore 
conceivable that the global field could result in an h/T 
scaling. 

All these considerations are only valid at the SU(2) 
invariant Bose-Fermi Kondo fixed point. However, 
magnetic interactions in CeCu5_ x Au x and most heavy 
fcrmion materials are strongly spin- anisotropic, and we 
have shown that anisotropy is relevant at the SU(2)- 
invariant Bose-Fermi Kondo fixed point. 

Thus if there is any locally critical quantum field the- 
ory that describes the QCP, it must be anisotropic. One 
of the candidates for the non-trivial fixed points is the 
new XY Bose-Fermi fixed point. At this fixed point, the 
coupling to the z-componcnt of the spin is irrelevant and 
scales to 0. This is therefore, a goad candidate to de- 
scribe systems like YbRh 2 Si2- x Ge x ,E3, where magnetic 
fluctuations are more of XY-type. 

In CeCu5_ x Au x , however, the magnetic degrees of free- 
dom are more Ising-like. The corresponding quantum 
phase transition can therefore possibly be described by 



the Ising Bose-Fermi fixed point of the local model with 
structure g x = g y = 0, g z ~ 1, and X z ^ \ x = X y . 



2. XY invariant theory 

As we discussed earlier, at this fixed point correlations 
in the xy and z spin directions decay with different pow- 
ers. This immediately implies that in a self-consistent 
theory one must also assume that the xy and z compo- 
nents of the field </>,; have also different anomalous di- 
mensions, e± and e z . It is straightforward to generalize 
our computations to this case. Similar to the e± = e z 
case, the asymptotic behavior of the correlations of the 
xy spin components is determined by a Ward identity, 
and is given by (S x (t)S x (0)) ~ (S v (t)S v (0)) ~ 
which yields again the exact relations, 



xy = 1 

e _L,QCP — 1 

7/ XY 

Vh _L 



1/2 



(19) 
(20) 



and implies conclusions similar to the ones drawn for the 
SU(2)-invariant fixed point for spin correlations within 
the xy plane. 

The other exponent, r] z , is however, not determined 
by a Ward identity And therefore the condition f] z (e z ) = 
2 — e z yields a non- trivial exponent e z ^ 1. In leading 
order we find that e z vanishes, and it next order it turns 
out to be negative, indicating that the leading terms 
probably originate from analytic corrections and scale 
as (S z (t)S z (0)) ~ 1/t 2 . These results are, however, only 
approximate, and not very reliable, since the e-expansion 
formulas are not applicable any more for ej_ = 1. 



3. Ising theory 

As we mentioned above, in CeCus_ x Au x the magnetic 
degrees of freedom are more Ising-like. The correspond- 
ing quantum phase transition is therefore most proba- 
bly be described by the Ising Bose-Fermi fixed point 
(dx = g y = 0, g z ^ 0, and X z ^ X x = X y ) within a 
dynamical mean field theory. This fixed point leads out 
of the region of applicability of our e expansion. Never- 
theless, from the Ward identity it still follows that 



UNI 
e z,QCP 

„UNI 

Vh,z 



I 



1/2, 



(21) 
(22) 



while the other exponents ej_ < 1 and r/± > 1 are now 
determined by the specific form of the beta functions. 

It is important to emphasize, that the prediction that 
for this locally critical dynamical mean field model local 
spin correlations in different directions decay with differ- 
ent powers and in one direction the decay is ~ 1/r while 
in some other Direction(s) it decays with a larger power, 
is rather robust, and it should be possible to test it by 
polarized neutron scattering. 
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II. BASIC DEFINITIONS 



B. Pseudofermion technique 



A. Spectral representations 

All calculations presented here were carried out in Mat- 
subara space using the spectral representation of the var- 
ious propagators. Though it is not necessary for the cal- 
culations, we found it rather useful and intuitive to con- 
struct also microscopic free Hamiltonians and fields that 
give rise to the local propagators in Eqs. (§) and (H|). 

We constructed the bosonic field 4> m terms of one- 
dimensional bosons, b a (k) as 

i- A 



dk\k\^^=(b a (k)+bi(-k)) 

A V 2 



(23) 



Here b a {k) destroys a boson with polarization a and mo- 
mentum k, and A is the high-energy cut-off. The 6 a (fc)'s 
are normalized such that [b a (k), bj 3 (k')} = 2ir 8 a p 5{k — 
k'). Using a free bosonic Hamiltonian to generate their 
dynamics: 



H 



(0) -E 

a=l 



r A dk 

J —\k\bt(k)b a (k), 



(24) 



it can be readily verified that Eqs. ( p3| ) and (p4j) lead 
to the propagator (||) with a constant T(2 — ey Note 
that cj) is not the usual free bosonic field corresponding 
to Eq. (]24|) that would be defined through 

<fP(x) ee [ A dk -J= (b a (k) + &t (-*)) . (25) 

It is also a trivial matter to show that the bosonic 
Matsubara Green's functions can be represented as: 

2? (0) (^)= (26) 



e(0 



(27) 



'-A 

sgn(0 Id 1 " 

with £>(£) the spectral function of the Bose field (f>. 

Similarly, to represent the fermionic propagator we can 
use free one-dimensional chiral fermions, described by the 
Hamiltonian 



H 



(0) 



J2 I VFk ^l<r^k,a 

„ — lJ-A z7r 



(28) 



with vf = 1 the Fermi velocity. The fermion field ip k a 
creates a fermion with momentum k, spin a, and energy 
k, and is normalized as {ip k CT , ipk'. a} = 2tt S(k— k'). Then 
the correlation function of the field 

(29) 



behaves exactly as Eq. ( 
be represented as 

G (0) (iw„) = 



dk , 

— Vk,a 



), and its Fourier transform can 



1 



_ A 2tt iuj n - £ 



(30) 



In order to treat the spin dynamics in a field the- 
oretical framework we used Abrikosov's pseudofermion 
technique E3 In this approach one represents the spin 
operator as S a -> J2 m ,m< fm S mm< fm> , where f m is a 
pseudofermion annihilation operator corresponding to 
the S z = m spin component. In terms of these pseud- 
ofermions, the interaction terms can be treated using 
standard field theoretic methods. However, in order to 
obtain a faithful representation of the spin one has to 
project out unphysical states with Y^mfmfm > 1- This 
is achieved by adding a term 



Hps — fio ^2 fmfn 



(31) 



to the Hamiltonian, taking the chemical potential /^o — > 
oo , and keeping only the leading contributions in the end 
of the calculation. 

The propagation of free pseudofermions is described by 
the Matsubara Green's function 

0<H(«*O = -(Tf m fl)t{ = . (32) 



III. RENORMALIZATION GROUP 

In order to analyze the behavior of the Bose-Fermi 
Kondo model at its non-trivial fixed points, we used field 
theoretical renormalization group methods, and applied 
expansion in the small parameter e. 

A. Spurious divergences and counterterms 

Simple power counting shows that the pseudofermion 
self-energy contains 'spurious' contributions that diverge 
as ~ A. These contributions only renormalize the chem- 
ical potential, /iq — ► /i, however, they lead to spurious 
divergences in higher order diagrams, and must there- 
fore be removed. 

To remove these divergences, we applied a counterterm 
procedure. In this approach we rewrite H ps as 



(33) 



and formulate the perturbation theory in terms of fi 
rather than /j,q: 



^ mm' 

iuj n - \i 



(34) 



The term ~ ^M^m fmfm must be treated as a pertur- 
bation and it systematically cancels all spurious diver- 
gences in higher order diagrams. The coefficient <5/i is 
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determined self-consistently up to a given order of the 
perturbation theory by the requirement: 

r< 2 ) = g- 1 (w n = n) = o , 

and is given in leading order by 
S(S + 1) A 



B. Vertex functions 



In order to compute the renormalization group equa- 
tions, in addition to the inverse pseudofermion propaga- 
tor T( 2 \ we have to compute the 'three-point' and 'four 



point' vertex functions, r« 



(3) 



and li 4) 



corresponding to 



the couplings 7 Q and A Q (see Fig. 0). This is straight- 
forward (though somewhat tedious) within the pseud- 
ofermion technique. 

As we shall see later, the fixed point couplings at the 
non-trivial fixed points turn out to be small, 7 Q ~ Q(\fe) 
and A Q ~ 0(e). This makes it possible to organize dia- 
grams according to the powers of e and perform a system- 
atic expansion. The leading order diagrams are shown in 
Fig.|. 

Having removed all 'spurious divergences' with the 
counterterm the diagrams contain only logarithmic in- 
frared singularities ~ In— in the e — * limit, and can 
be expanded as a power of the logarithm. In leading 
order, the vertex functions are given by the following ex- 
pressions: 



r< 2 >=u>{l + lnA [ij>] + ...} 

1 1 a 

+ ln R [\(9x ~ 9y ~ 9z)] +■■•} 

4 4) =K(i-- A {g x -9y-g z )) 



. A 
+ ln— 

M 



(35) 



(36) 



(37) 



To perform the computation up to second order in e 35 
diagrams must be taken into account. Their list and 
the corresponding rather lengthy results are given in the 
Appendix. 



C. The renormalization group transformation 

As mentioned above, once we removed all spurious di- 
vergences, the perturbation series contains only logarith- 
mic divergences in the e — > limit. These logarithmic di- 
vergences can be systematically summed up and treated 
using renormalization group methodsfia 



The chemical potential fi has, of course, no physical 
meaning. Therefore, once the unphysical contributions 
corresponding tojpultiple occupancy of the /-states have 
been eliminated,E3 we can transform it out by the trans- 
formation iu> n — ► u> + fi. After this transformation, all 
quantities depend exclusively on A, the dimensionless 
couplings, and the ratio uj/A. 

To proceed further we define the pseudofermion Z- 
factor as 



duj 



Z-^u/A, {%}] 



(38) 



where {qt} (i = 1,..,6) refers to the six dimensionless 
couplings of the theory, and l^ 2 ) denotes the inverse 
pseudofermion propagator. 

Then the renormalizability of the theory implies that 
there exists a transformation A — ► A', — > q^ such that 
the 'renormalized' propagators and vertices remain in- 
variant for all frequencies up to G(w/A) corrections: 



r 



(«) - 



Z(w/A; ft ) T^{uj,A- qi ) = fixed 



(39) 



where T^ 2 -* has been defined earlier, and n = 3 and n — 
4 refer to the bosonic and fermionic vertex functions, 
respectively (see Fig. ||). 

This invariance is expressed through the Callan- 
Symanzik equation: 



dl 



-inr 



d_ 

dl' 



inr 



(n) 



E 



dgj d 
dl dqi 



inr 



(n) 



, (40) 



where we introduced the scaling parameter I = ln(Ao/A), 
with Ao the bare value of the cut-off. 

Substituting Eq. @ into Eq. @ we obtain the fol- 
lowing differential equations: 



<91nr( J 
dl 



d\nZ 
dl 



^ dq a /d\nT^ 
^~aT \ dq a 



dlnZ\ 
dq a J 



(41) 



To derive the scaling equations, we take the e — ► 
limit. Note that Eq. (|4l]) is satisfied for all frequencies 
and temperatures provided that they are much smaller 
than the cut-off A. Therefore, expanding all vertex func- 
tions in powers of ln(-^j), the Callan-Symanzik equations 
can be obtained by just picking the constant terms and 
the terms proportional to ln(^j). [Of course, in a renor- 
malizable theory the same equations must be generated 
by comparing higher order terms in th expansion. In- 
deed, as a check we verified in leading order in e that 
the same scaling equations are obtained from comparing 
the ~ ln(^j) and ~ In (^) terms too.] Plugging in the 
expressions of Eqs. in the Eq. (0) we obtain the 

scaling equations Eq. and Eq. (|g). The second order 
C(e 2 ) scaling equations are given in the Appendix. 
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<(3) 



+ 



derivative 5Y,/5Q, and thus the bosonic vertex function 
by the Ward identity: 



r 



(2) 



.— = a- £/2 — r (3) 

dh a -i a a 



(46) 



Another way to derive this equation is to notice that 
g(o) 

can be expanded as 



,(4) 



+ 




+ 



FIG. 4. Leading order corrections for the pseudofermion 
self-energy and vertex functions. Dashed lines stand for 
free pseudofermion propa gat ors [Eq. (^)] wavy lines denote 
bosonic propagators [Eq. (|26|)1 and continuous lines represent 
fermion propagators [Eq. Q3C|)] . 



D. Local magnetic field 

1. Scaling equations 

To investigate the effect of local magnetic field we have 
to add the term Eq. (^) to the Hamiltonian. As a conse- 
quence, £/(°) is modified to: 



7(0) _ 



(42) 



and S also becomes a matrix in spin space. Restricting 
the calculation to the case where h represents the small- 
est energy scale in the problem, \h\ < uj,T, we can ex- 
pand the renormalized inverse pseudofermion propagator 



as r 



(2) 



ee z [r : 



(2) 



KT^l + 0(h 2 )} , 



(43) 



(2) 

where r o and Z denote the inverse propagator and Z- 

(2) 

factor in the absence of magnetic field, and T\ ' a is also a 
matrix in spin space: 



-.(2) 



8T, 
dh a 



h=0 



(44) 



From the invariance of T R it immediately follows that 
the magnetic field satisfies the following scaling equa- 
tions: 



7(0) 



IUJ — /i 



S a - + ... (47) 

10J — jJL IU> — jl 



Differentiating £ with respect to h a at h = thus 
amounts in diagrams where one of the pseudofermion 
lines has been replaced by an external vertex j a (see 
Fig. ||). It is easy to see that each vertex diagram is 
generated in this way exactly once and with the cor- 
rect weight. Combining the Ward identity Eq. (jitj ) 
with the invariance of the renormalized vertex function, 
4j\^{Z r' 3 )) = 0, we immediately obtain that 



dl 



t m 



VA e / 2 7 J 







(48) 



which, taking into account the definition of the (3 function 
leads to the scaling equation: 



1 dh a 
h„ dl 



(49) 



or equivalently, Eq. (^) for the dimensionless magnetic 
field, h = h/A. 



4> 



i 



i 



FIG. 5. Diagrammatic representation of the Ward iden- 
tity, Eq. (ph. 



E. Susceptibility 

There are two possibilities at a non-trivial fixed point 
with couplings g a = g* a an d A Q = A* : 
(a) If g* a f then (g* , A*) = and from Eq. (g) it 
follows that 



dl 



In 



Zh a Y\ 



(2) 



= . 



(45) 



However, E only depends on h through the propagator 
Q(°\ Therefore T^ a is proportional to the functional 



dh Q 
~df 



x = (1-J)fc„ 



(50) 



This immediately implies the exact result for thei-a com- 
ponent of the spin correlation function at T = OEJ 
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(TS a (T)S a (0)) 9 *^ = 



(r* T y ' 



if 5 *^0, (51) 



with T* a dynamically generated energy scale that is in- 
variant undep,the scaling transformation, similar to the 
Kondo scale.cJ 

On the other hand, if 

g* = then Hm g .— {/% ■ (g*,A*)/g*} 7^ takes some 
non-trivial value at the fixed point and no such universal 
relationship exists, and the spin-spin correlation function 
decays as 



(TS a (r)S a (0)) g ,^ 



with an exponent 



1 



if <£ = (), (52) 



(53) 



The spin correlation functions decay as a power law 
only at T — temperature. It is, however, easy to ex- 
tend these results to finite temperature by assuming that 
there exists a boundary conformal field theoryEj that cor- 
responds to the critical dynamics of the fixed point. In 
this boundary conformal field theory the impurity lives 
at the x = line of the complex plane, z = r + ix, and 
the theory is invariant under conformal mappings that 
map this boundary onto itself. 

The finite temperature correlation functions can be ob- 
tained by mapping the complex plane on a strip of width 
[3 = 1/T using e.g. the function w = ^artg z. Assum- 
ing that S a transforms as a primary field of dimension 
W2 



(S a (w 1 )S a (w 2 )) T = 

V dw\ ) V dw2 



(S a ( Zl )S a (z 2 )) , (54) 



we obtain that at the critical point: 

ttT 



Xa(r,T)= ( 



T*sin(7rTr) 



(55) 



Taking then the Fourier transform of Eq. ( p5[ ) and con- 
tinuing it analytically to the real axis one obtains an es^ 
act expression of the finite temperature susceptibility:E3 



T e it*- 1 / 2 r(i 



t r(t)r(|-t 



2ttT 



(56) 



The scaling form Eq. (|13J) readily follows from the asymp- 
totic properties of the T function. 



IV. CONCLUSIONS 



found that in agreement with Ref. [12], spin anisotropy 
is relevant at the purely bosonic and Bose-Fermi fixed 
points. 

We found two types of new bosonic fixed points: One 
one of them has XY symmetry, and is unstable against 
breaking the XY symmetry This fixed point may be 
relevant for non-magnetic impurities in antiferromag- 
nets at their quantum critical point with easy plain 
anisotropy, and it may also describe the coupling of the 
spins to the Goldstone modes in an ordered phase of an 
antiferromagnetES. The other fixed point has an Ising- 
like structure, and is stable. This fixed point could be 
relevant for non-magnetic impurities in an Ising antifcr- 
romagnet at the quantum critical point. 

We identified furthermore two new Bose-Fermi quan- 
tum critical points in the impurity model with XY and 
Ising symmetries, respectively, which parallel the bosonic 
fixed points and separate the bosonic and fcrmionic quan- 
tum phases. We have shown that local spin correlations 
at these fixed points decay as {TS Xj y(T)S Xi y(0)) ~ 1/t £ 
and (TS z (t)S z (0)) - 1/t £ at the XY and Ising fixed 
points, respectively, as guaranteed by a Ward identity. 
Correlations in the other directions decay faster with 
non-universal exponents that we computed to second or- 
der in e. These fixed points may be relevant at the 
locally quantum critical p puit s of YbRh2Si2- x Ge x and 
CeCu5_ x Au x , respectively.ll3o 

As discussed in detail in the Introduction, these re- 
sults have very important implications for a possible dy- 
namical mean field theory of the locally quantum critical 
behaviorLij: 

(i) First of all, they imply that the real part of the lo- 
cal susceptibility diverges logarithmically at the quan- 
tum critical point in the easy axis/easy plain directions, 
Xoasy, local H ~ m ( T *M> wh i\c the imaginary part is 

singular, XcasyjocaiM ~ s S n M- In tnc otner directions, 
however, the local susceptibility should show a power-law 

behavior, xLot- easy, local M ~ Xaot-easy.localH ~ 

ujP , and is likely to be dominated by analytic contribu- 
tions with a = and 0=1. These predictions should be 
directly observable by polarized neutron scattering and 
NMR relaxation measurements. 

(ii) A logarithmically divergent local susceptibility is con- 
sistent with the experimental scaling form, Eq. jT^), 
only if the magnetic fluctuations are effectively two- 
dimensional. This quantum critical behavior can there- 
fore be only approximate, and if it applies to the systems 
above, then a cross-over to some other type of quantum 
critical behavior should occur at lower temperatures. 

(iii) Finally, the above results imply an h 2 /T scaling for 
the local magnetic field. Whether this is consistent with 
the experimentally observed h/T scaling with the global 
field, is an open question that needs to be answered. 



In the present paper we analyzed the anisotropic Bose- 
Fermi Kondo problem using e-expansion methods. We 
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V. APPENDIX 




FIG. 7. Second order corrections in e for the bosonic ver- 
tex function. 



In this appendix we give a few detailed expressions that 
may be important to reproduce or extend our results. 

The easiest quantity to compute is the pseudofermion 
self-energy, whose diagrams are given in Fig | to C(e 2 ), 
and correspond to the following expression: 

r(2) « «{l + (J + ~ E InA} , 

(57) 

Note that the counterterm itself is also proportional to e 
and is therefore important only in the 0(e 2 ) calculation. 




FIG. 6. Second order corrections in e for the pseud- 
ofermion self-energy. The cross denotes the counterterm. 

The expression of the bosonic vertex is somewhat 
lengthier: 

ri 3) = A £ / 2 7Q {^o Q (ffa, A Q ) + \n^-F?(g a ,\ a ) + ...}, 

(58) 

where the coefficients Fq and F\ are given by 

F% = 1- \(g x -9y-g z ) + 0(e 2 ) , (59) 

Fx = ^~^(9x ~9y- 9z) ~Jq(91 + 9 2 y + 9z) 

3 , .7 
+ ^9x{9z+9v) - g9y9z 

+ \{\l-\ 2 y -\l) + 0{e 2 ). (60) 
The corresponding diagrams are depicted in Fig. 



Finally, the vertex corrections to the usual Kondo in- 
teraction to this order are depicted in Fig. |j| and the 
corresponding expression reads: 

lt 4) = Q?(ffa, A a ) + \n—Q1{g a , A a ) + . . . , (61) 

where Qo and Qi are given by 

Q x = K(i-\{g x ~g y -g z ) + O(e 2 )) , (62) 

Of = XyX z (1 + \{g y + g z - g x )) + K Ff . (63) 

Here we only gave the x-component of the vertex func- 
tions, the others can be obtained by cyclic permutation. 




FIG. 8. Second order corrections in e for the fermionic 
vertex function. 



The scaling equations can then be simply obtained by 
substituting these expressions into Eq. (|4l]) and are given 
by 

= \K - ^\ x (g y + g z ) - ^(A 2 + A 2 ) 
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*x(gx(g y + gz) - g y gz) + C(e 4 ) , 
9z- g v ) + gl(g y + 9z) 



Fixed point 


{Aa} 


{3<*} 


A 


= -Hi 




symm 

yr 


R iulp (T) 


Stability 


Kondo SU(2) 


X a — A — > oo 


3a = 




1 


none 


none 


Ci - c* 2 r^ 


STABLE 


SU{2)Bose 


A a = A / 


g a = o 


2 




e , e 2 
- 4- — 

2^2 


stable 




unstable 


5(7(2) Bose-Fermi 


X a = A / 


Sc. = # 


e 


12 


1 + 2£ 

2 ^ 24 


2 ' 6 


Ci + C 2 T A 


unstable 


XY Bose 


A Q =0 


9x = 3a ^ 3z = 




6/2 


6 + 6 2 


stable 


/ji2A 


unstable 


XY Bose-Fermi 


A x = Ay 7^ A z 


3z = 3y + 9* = 


e 


9e 2 
8 


fc T 4 


0.618e + 0.123e 2 


Ci + c 2 r A 


unstable 


Uniax Bose-Fermi 


A x = A y 7^ A z ~ 1 


3a: = 9y = 0, g z ~ 1 




? 


? 


? 


? 


unstable 


Uniax Bose 


A Q =0 


3* = 3a - 0, 3^ -» 00 




? 


none 


stable 


? 


STABLE 



TABLE I. Stability and structure of the various fixed points in the absence of magnetic field: (a) y,: is the dimension of 
the leading irrelevant operator. This exponent governs thermodynamics and scattering rates at the fixed point, (b) y r is the 
dimension of the leading relevant operator. This corresponds to breaking the symmetry of the fixed point, (c) y* ymm is the 
scaling dimension of the leading relevant operator within the subspace of given symmetry. This governs the quantum phase 
transition between the Bosonic and Fermionic Kondo fixed points. 
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